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Abstract. Given a finite dimensional manifold A^, the group Diff5(Af) of 
difleomorphism of A^ which fall suitably rapidly to the identity, acts on the 
manifold B(M, N) of submanifolds on N of diflfeomorphism type M where M 
is a compact manifold with dim A/ < dimAf. For a right invariant weak 
Riemannian metric on Diff5(Af) induced by a quite general operator L : 
^s{N) ~^ r(T* Acgivol(A)), we consider the induced weak Riemannian metric 
on B{M, N) and we compute its geodesies and sectional curvature. For that 
we derive a covariant formula for curvature in finite and infinite dimensions, 
we show how it makes O'Neill's formula very transparent, and we use it finally 
to compute sectional curvature on B{M,N). 



1. Introduction 

It was 46 years ago that Arnold discovered an amazing link between Euler's 
equation for incompressible non-viscous fluid flow and geodesies in the group of 
volume preserving diffeomorphisms SDiff(IR") under the L^-metric [2]. One goal 
in this paper is to extend his ideas to a large class of Riemannian metrics on the 
group of all diffeomorphisms Diff5(iV) falling suitably to the identity, of any finite 
dimensional manifold N . The resulting geodesic equations are integro-differential 
equations for fluid-like flows on N determined by an initial velocity field. In previous 
papers [121 HOl H] : we have looked at the special case where N = R" and the metric 
is a sum of Sobolev norms on each component of the tangent vector but here we 
develop the formalism to work in a very general setting. 

The extra regularity given by using higher order norms means that these metrics 
on the group of diffeomorphisms can induce a metric on many quotient spaces of 
the diffcomorphism group modulo a subgroup. This paper focuses on the space 
of submanifolds of N diffeomorphic to some M which we denote by B{M,N). 
Diff5(A^) acts on B{M, N) with open orbits, one for each isotopy type of embedding 
of M in A''. The spaces B may be called the Chow manifolds of N by analogy with 
the Chow varieties of algebraic geometry, or non-linear Grassmannians because 
of their analogy with the Grassmannian of linear subspaces of a projective space. 
The key point is that the metrics we study will descend to the spaces B{M, N) so 
that the map Diff5(iV) B{M,N) (given by the group action on a base point) 
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is a Riemannian submersion. Geodesies from one submanifold to another may be 
thought of as deformations of one into the other reahzed by a flow on N of minimal 
energy. 

In the special case where M is a finite set of points, B{M,N) is called the 
space of landmark point sets in N. This has been used extensively by statisticians 
for example and is the subject of our previous paper [8]. The case B{S^,M.'^) is 
the space of all simple closed plane curves and has been studied in many metrics, 
see [71[TT1[T3] for example. This and the case _B(S'^,R'^) of spheres in 3-space have 
had many applications to medical imaging, constructing optimal warps of various 
body parts or sections of body parts from one medical scan to another fl3[ 117] . 

The high point of Arnold's analysis was his determination of the sectional curva- 
tures in the group of volume preserving diffeomorphisms. This has had considerable 
impact on the analysis of the stability and instability of incompressible fluid flow. 
A similar formula for sectional curvature of B{M, N) may be expected to shed light 
on how stable or unstable geodesies are in this space, e.g. whether they are unique 
and effective for medical applications. 

Computing this curvature required a new formula. In general, the induced inner 
product on the cotangent space of a submersive quotient is much more amenable to 
calculations than the inner product on the tangent space. The first author found 
a new formula for the curvature tensor of a Riemannian manifold which uses only 
derivatives of the former, the dual metric tensor. This result, Mario's formula, is 
proven in section 2. In this section we also define a new class of infinite dimensional 
Riemannian manifolds, robust Riemannian manifolds to which Mario's formula and 
our analysis of submersive quotients applies. We also obtain a transparent new 
proof of O'Neill's curvature fomula. This class of manifolds builds on the theory 
of convenient infinite dimensional manifolds, see [6]. To facilitate readability this 
theory is summarized in an Appendix. 

In section 3, we describe diffeomorphism groups of a finite dimensional manifold 
N consisting of diffeomorphisms which decrease suitably rapidly to the identity on 
A'^ if we move to infinity on N; only these admit charts and are a regular Lie groups. 
We shall denote by Diff 5 ( A^) any of these groups in order to simplify notation, and 
by Xs{N) the corresponding Lie algebra of suitably decreasing vector fields on A^. 
We introduce a very general class of Riemannian metrics given by a positive definite 
self-adjoint differential operator L from the space of smooth vector fields on A^ to 
the space of measure- valued 1-forms. This defines an inner product on vector fields 



Note that LX paired with Y gives a measure on A^ hence can be integrated without 
assuming A^ carries any further structure. Under suitable assumptions, the inverse 
of L is given by a kernel K{x,y) on N x N with values in p\TN ®P2TN. We 
then describe the geodesic equation in Diff5(A^) for these metrics. It is especially 
simple written in terms of the momentum. If ip{t) £ Diff 5 (A?^) is the geodesic, then 
X[t) = 9t(<p) o ^p~^ is a time varying vector field on A^ and its momentum is simply 



X,y by: 




LX{t). 
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In section 4 we introduce the induced metrics on B{M, N). We give the geodesic 
equation for these metrics also using momentum. One of the keys to working in 
this space is to define a convenient set of vector fields and forms on B in terms of 
auxiliary forms and vector fields on N. In this way, differential geometry on B can 
be reduced to calculations on N. Lie derivatives on N are especially useful here. 

In the final section 5, we compute the sectional curvatures of B{M,N). Like 
Arnold's formula, we get a formula with several terms each of which seems to play 
a different role. The first involves the second derivatives of K and the others are 
expressed in terms which we call force and stress. Force is the bilinear version 
of the acceleration term in the geodesic equation and stress is a derivative of one 
vector field with respect to the other, half of a Lie bracket, defined in what are 
essentially local coordinates. For the landmark space case, we proved this formula 
in our previous paper [8| . We hope that the terms in this formula will be elucidated 
by further study and analysis of specific cases. 



2. A CovARiANT Formula for Curvature 

2.1. Covariant derivative. Let {M,g) be a (finite dimensional) Riemannian 
manifold. There will be some formulae which are valid for infinite dimensional 
manifolds and we will introduce definitions for these below. For each x € M we 
view the metric also as a bijective mapping ■ T^M T*M. Then g^^ is the 
metric on the cotangent bundle as well as the morphism T*M TM . For a 1-form 
a e n^{M) = r{T*M) we consider the 'sharp' vector field a» = g^^a e X{M). If 
a = Uidx^ , then a" = atg^-^dj is just the vector field obtained from a by 'raising 
indices'. Similarly, for a vector field X G X{M) we consider the 'fiat' 1-form 
= gX. U X ^ X'di, then X^ = X'gijdx^ is the 1-form obtained from X by 
'lowering indices'. Note that 

(1) a(/3«) = g-'{a,P)^g{aKl3^) = /3(a«). 

Our aim is to express the sectional curvature of g in terms of a, j3 alone. It is 
important that the exterior derivative satisfies: 

(2) da(/3«,7») = (/3*)a(7») - (7")a(/S») -a([/3",7»]) 

We recall the definition of the Levi-Civita covariant derivative V and its basic 
properties: 

(3) 2g{VxY,Z) 

(4) (Vxa)(y) 
g{{Vxa)\Y) 

(5) V;f(a«) 
Xg-\a,P) 



= X{g{Y, Z)) + y(.g(Z, X)) - Z{g{X, Y)) 

- g{X, [r, Z]) + g{Y, [Z, X]) + .g(Z, [X, Y]) 
= Xa{Y) - a{VxY) 
= {Vxa){Y) = X{a{Y)) - a{V xY) 
= Xg{a\Y) - g{a\VxY) = g{Vxa\Y) 
= (Vxa)» 

= g-\Vx(y,fi)+g-\a,SIxP) 
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From this follows 

2(V„./?)(7») = 25-i(V„.^,7) = 25((V„,/3)»,7*) = 2ff(V„,/3«, 7") = 

= ip, 7) + P^g-' (7, a) - 7«5"' (a, 

- 9-i(a, [^«,7«]') + 9-\P, [7«,a«]^) +5-^(7, 

= Q!*'/3(7») + /3*'7(af ) - 7*'^(a'*) 

-a([/3«,7"])+/3([7",a«])+7([a«,^«]) 

= /3»7(a*') - a([/3*', 7I']) + j{[ai, p^) - dj3{-fK a") 

(7) = a«7(/3«) - /3*'a(7*') + 7*'a(/3*') + da(^'',7*') - d/3(7*',a*') - d-f{a\P^) 

2.2. Theorem. (Mario's Formula) Assume that all 1-forms a, /3, 7, 5 G ^^(M) 
are closed. Then curvature is given by: 



g(i?(a«,/?«)7«,<5«) =i?i+i?2 + i?3 

R2 = i (-<7-Hrf(7(/3«)),d(<5(a«))) +.9-i(d(7(«»)), 

i?3 = i (.g([<5«,««],[,5»,7^]) -.9([<5«.„5^],[a«,7*]) +25([a«,/3«],[7",'5«])) 



For the numerator of sectional curvature we get 



g{R{J,(3^)(3^,J) i?i + i?2 + i?3 

i?i = i (a«a«(||/3f ) - + /?«a«).9-i(a,/3) +/3«/3»(||Qf)) 

^2 = i (|M(ff-i(a,/3))f - g~\d{\\af),d{\\f3f)) 
fi3 = -|-5(M,/3"],M,/3"]) 



Recall that sectional curvature is then 

g(i?(a«,/3«)/3«,a») 



A;(a«,/3«) 



^||/3P-5-i(a,;3)2 



Proof. We shall need that for a function / we have: 
(V;3«7)"/ = dfiiV^n)^) = g-\df,y0n) = Ph-\df,^) - g-\Vp,df,^) 

= Ph^f - V0«#(7*) = Ph^f - 5/3"#(7") + i#"7(/5") - h^m^) 
(8) = id/«7(/3«) + ^[/5«,7"]/ = id(7(/3»))(rf/«) + ^[/3«,7«]/ 

For the three summands in the curvature formula, by multiple uses of formulas (2) 
and (7) and the closedness of a, /3, 7, 5, a straightforward calculation gives us: 

4(V„,V^.7)(^") = 

= 2a«(V^,7)('^") - 2(V^,7)"'5(a«) + 25«a((V^,7)«) - 2d(V^,7)((5«, a«) 
= 2a«(V;3,7)(<5") - d{-i{P^)){d{5{a^f) - [/3», 7»](5(a») 

+ 2a«(V^«7)(^") + 2(V;3«7)([^", a']) = ■ • ■ = 
= -g-^ (d(7(/3«)), rfl-^la"))) + 
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and similarly 

- 4(V^. V„.7)(<5«) = +g-\d{^{a% d{6{pi))) - g{[SK 7"]) 

- 2l3^a^'y{5^) + 2P^-f^S{J) - P^S^ai-y^) + [a", 7*']5(/3*') - (5*'/3*'a(7*') 
-2(V[,«,^,]7)('5") = 

= -[a", /3*']7(^'') + T^^Cla", ^f]) - 5h{[aK Z?"]) - rf^, ^'']^(7f , 5") 

Now we can compute the curvature (remember that da — dP = ■ ■ ■ = 0): 

45(i?(a«,/3«)7«,<5«) = 4J(ii(a«,/3«)7«) = 4^(V„,V^«7« - V^jV^a" - V[„»,^»]7*) 
= 4(V„, V^,7 - V^, V,«7 - V[„t,,,t]7)((5«) 
= -g-\d{^{(3i)),d{6{J)))+g-'{d{j{Jj),d{S{P% 

+ g{[Si,a% 7«]) - g{[S\l3^], 7"]) + 25([a«, [7", <5«]) 

- aS^-^l/?") + a*(5*/3(7*) + /3''7*<5(a«) - ;3'',5''a(7'') 

For the sectional curvature expression this simplifies (as always, for closed 1-forms) 
to the expression in the theorem. □ 



2.3. Mario's formula in coordinates. The formula for sectional curvature 
becomes especially transparent if we expand it in coordinates. Assume that a = 
aidx^ , /3 = Pidx^ where the coefficients ai,f3i are constants, hence a,/3 are closed. 
Then a" = g^^aidj,(3'^ = g^^(3idj. Substituting these in the terms of the right hand 
side of Mario's formula for sectional curvature, we get: 

2nd deriv. terms = 2i?i = 2JaH}\l3f ) + 2/3^/3H\\af ) - 2{a^/3^ + P^J)g-\a,/3) 

= 2a,g'''{ajg^'{Mg''').t),s + 2/3.5'^(/3j.g^'*(afea/5"),*),. 

- 2aig*^(/3,g^*(/3feaz/'),0,s - 2Ag^^(a,.9^*(«,A/'),t),a 

= 2{aA - Oikl^^) ■ {ajpi ~ aip,) ■ 5*^(5^*5^/),. 
1st deriv. terms = 4ii2 = ||%-'(a,^))f-ff-im|af),rf(||/3f)) 

= (a,/3j5'^),,5'*(a,/3fc5''),t - {aiajg'^),sg'\Mig'''),t 

= -\{aiPk - akPi) ■ {oijPi - aipj) ■ g^g'^g'^i 
Lie bracket = [a«,/3»] = [aig'^ {^kg^') ,s - Pig'' {akg^') ,s) dt 

= {aipk - akl3i)g''g'^:dt 
Lie bracket term = 4R3 = -3g{[aKp% [a",^"]) 

= -3{aA - akPi) ■ {ajpi - ai/3j) ■ g'' g*;! g^^g^' g\l 
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hence we have the coordinate version for the three terms in sectional curvature: 



Note that the usual contravariant metric tensor gij occurs in only one place, ev- 
erything else being derived from the covariant metric tensor g''^ . Note that the fist 
term Ri can be split into a pure second derivative term Ru = g^^g^^g'^lt plus a first 
derivative term R12 = g^^g^^lg^i- 

There is also a version of Mario's formula which is, in a sense, intermediate 
between the coordinate free version and the coordinate version. The main thing 
that coordinates allow you to do is to take derivatives using the associated flat 
connection. In the case of this formula, this introduces auxiliary vector fields Xa 
and Xp playing the role of 'locally constant' extensions of the value of a" and /3' at 
the point x G M where the curvature is being calculated and for which the 1-forms 
a, /3 appear locally constant too. More precisely, assume we are given Xa and X^ 
such that: 

(1) X„(X)=Q«(X), X^{x)=PKx), 

(2) Then a" — X^ is zero at x hence has a well defined derivative Dx{a} — X^) 
lying in Hom(Tj,M, T^M). For a vector field Y we have D.^{a^ - Xa).Yx = 
[Y, a« - Xa]{x) = Cvia* - Xa)\x- The same holds for /3. 

(3) Cx^ (a) = Cx^ (/3) = Cx, (a) = Cx, (/3) = 0, 

(4) [X^,Xp]=0. 

Locally constant 1-forms and vector fields satisfy these properties. Using these 
forms and vector fields, wc then define: 

F(a,/3) : = ^d{g~^{a, P)), a 1-form on M called the force, 

D{a,p){x) : = D4l3^ - Xp).aHx) 

= d{(3^ — Xi3).a^(x), a tangent vector at x called the stress. 

Then in the notation above: 

(i?(a«, a») (x) = Rn + R12 + R2 + R3 

^11 = 5 (^1 i9-'){P, P) - '^Cx^Cx, {g-'Ka, 13) + (.9-i)(a, a)) 

i?i2 = (i^(a,a),D(/3,/3)} + (F(/3, /?), i?(a, a)) - (F(a, D(a, /3) + D(/3, a)) 

i?2 = - (F(a,a)),F(;3,;3))^_,) {x) 

fi3 = -fP(a,/3)-mQ)llL 

The reformulation of R\ follows from the calculation: 

aW(||/3||2)(a;)=X„a»(ll^f)(^) 

= X^XM?){^) + - ^a)(Pf )(^) 



SOBOLEV CURVATURE 



7 



= Cx^CxA9'\P,P)){x) + {Dx{a^ - X^).X^{x),d\\pf ){x)) 
= (5-i)(/3, /3)(a:) + - X^).XAx), d\W){x)) 

and the similar result for the other terms. The reformulation of i?3 comes from the 
calculation: 

[a»,/3»](x) = (X„ o /3«)(x) - (X^ o J){x) 

= {X^ o (/3« - X^))(a;) - {Xfi o (a« - 
= - X^).X„(a;) - - 

2.4. Infinite dimensional manifolds. The main focus of this paper are the 
infinite dimensional manifolds of diffeomorphisms of a finite dimensional N ^ of the 
embeddings of one finite dimensional M into another X and of the set of subman- 
ifolds of a manifold X . These are infinite dimensional and can be realized in 
multiple ways depending on the degree of smoothness imposed on the difFeomor- 
phism/embedding/submanifold. The first two have realizations as Hilbert man- 
ifolds but the last does not. Moreover, the group law on the Hilbert manifold 
version of the group of diffeomorphisms is not differentiable. If one desires to carry 
over finite dimensonal techniques to the infinite dimensional setting, it works much 
more smoothly to use the Frechet space of C°° functions decreasing rapidly at in- 
finity as the base vector space for charts of these spaces. But then its dual is not 
Frechet, so one needs a bigger category for charts on bundles. The best setting has 
been developed by one of the authors and his collaborators [(3, and uses 'c°°-open' 
subsets in arbitrary 'convenient' locally convex topological vector spaces for charts. 
This theory and some of the reasons why it works are summarized in the appen- 
dix. For our purposes, complete locally convex topological vector spaces (which are 
always convenient) suffice and, on them 'c°°-open' just means open. 

To extend Mario's formula to infinite-dimensional manifolds then, let (M, g) be 
a so-called 'weak Riemannian manifold' a convenient manifold M and smooth 
map: 

g : TM Xm TM — > M 

which is a positive definite symmetric bilinear form on each tangent space 
TxM^x e M. For a convenient manifold we have to choose what we mean by 
1-forms carefully. For each x e M the metric defines a mapping g^ : T^M — > T*M 
(which we denote by the same symbol gx)- In the case of a Riemannian Hilbert 
manifold, this is bijective and has an inverse but otherwise is only injective, hence 
the term 'weak metric'. The image g{TM) C T*M is called the g-smooth cotangent 
bundle. Then g^^ is the metric on the f/-smooth cotangent bundle as well as the 
morphism g{TM) TM. Now define rjJ(M) := T{g{TM)) and a" = g'^a e 
X{M)^X^ = gX are as above. The exterior derivative is now defined by: 

da{p\-i^) = (/3")a(7«) - (7«)a(/3«) - a([/3'', 7"]) 

We have d : rji(M) ^ n'^{M) = r(L2j^^„(TM; M)) since the embedding g(TM) C 
T*M is a smooth fiber linear mapping. Note that on an infinite dimensional mani- 
fold M there are many choices of differential forms but only one of them is suitable 
for analysis on manifolds. These are discussed in [H Section 33]. Here we consider 
subspaces of these differential forms. 
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Further requirements need to be imposed on (Af , g) for our theory to work. Since 
it is an infinite dimensional weak Riemannian manifold the Levi-Civita covariant 
derivative might not exist in TM . Existence of the Levi-Civita covariant derivative 
is equivalent to the condition that the metric itself admits gradients with respect 
to itself in the following senses. The easiest way to express this is locally in a chart 
U C M. Let Vij be the vector space of constant vector fields on U. Then we assume 
that there are smooth maps gradig and grad2(7 from U x Vu >^ Vu to Vjj, bilinear 
in Vu X Vu such that 

D.,zgx{X, Y) = g,{Z, gradi gix){X, Y)) = g^{gv&d^ g{x){Z, X), Y) 

(If we express this globally we also get derivatives of the vector fields Y and Z.) This 
allows to use (12.11 3) to get the covariant derivative. Then the rest of the derivation 
of Mario's formula goes through and the final formula for curvature holds in both 
the finite and infinite dimensional cases. 

Some constructions to be done shortly encounter a second problem: they lead to 
vector fields whose values do not lie in T^M, but in the Hilbert space completion 
TxM with respect to the inner product g^- To manipulate these as in the finite 
dimensional case, we need to know that [J^^pj T^M forms a smooth vector bundle 
over M. More precisely, choose an atlas {Ua, Ua : Ua ~^ E) of M, where the Ua C 
M form an open cover of M, where each Ua '. Ua ^ E is a, homcomorphism of Ua 
onto the open subset Ua{Ua) of the convenient vector space E which models A/, and 
where Ua0 = Ua° '■ Ufi{Ua n Up) — >■ Ua{Ua n C/^)) is a smooth diffeomorphism. 
The mappings x n- (pap{x) — duai3{u'^^{x)) G L{E,E) then form the cocycle of 
transition functions (pai3 Ua H Up ^ GL{E) wich define the tangent bundle TM . 
We then assume that the local expression of each Riemannian metric gx on E 
are equivalent weak inner products hence define Hilbert space completions which 
are quasi-isometric via extensions of the embeddings of E (in each chart). Let 
us call one such Hilbert space %. We then require that all transition functions 
'^ap{x) : E ^ E extend to bounded linear isomorphisms % ^ % and that each 
'■Pap Ua ri Uj3 ^ L{H, W) is again smooth. 

These two properties will be sufficient for all the constructions we need so we 
make them into a definition: 

Definition. A convenient weak Riemannian manifold (M, g) will be called a robust 
Riemannian manifold if 

(1) The metric gx admits gradients in the above two senses, 

(2) The completions T^M form a vector bundle as described above. 

Note that a Hilbert manifold is automatically robust. We can make the relation- 
ship between robust manifolds and Hilbert manifolds more explicit if we introduce 
another definition, that of a pre Hilbert manifold similar to the notion of a pre- 
Hilbert topological vector space: 

Definition. A robust Riemannian manifold (Af, g) is called pre-Hilbert if there 
exists an atlas (Ua,Ua : Ua ^ E) for which: 
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(1) Each Ua{Ua) is contained in the Hilbert norm interior of its closure in H, 
which we denote Ua{Ua)^- 

(2) Ah chart change maps Mq^ extend to smooth mappings between the open 
subsets Ua{Ua)^, hence define a completion M C which is a Hilbert 
manifold. 

Note that in this definition the atlas must be properly chosen: for example its 
open sets Ua must be open in the weak topology defined by path lengths. More 
precisely, for any weak Riemannian manifold AI, the inner products assign a 
length to every smooth path in M and we get a distance function d{x, y) as the 
infimum of lengths of paths joining x and y (which might however be zero for some 
X ^ y). The topology defined by path lengths is usually much weaker than the 
strong topology given by the definition of M. 

These distinctions are well illustrated by the spaces we will discuss below. Firstly, 
manifolds of smooth mappings like Emb(M, N) with their canonically induced 
Sobolev metrics of order s > dim M/2 do admit completions Emb'*(M, N) to Hilbert 
manifolds hence are pre-Hilbert; see [61 42.1] for the explicit chart changes. But 
their quotient manifolds B{M,N) — Emb(M, iV)/ Diff(Af) are only robust in gen- 
eral because the second condition fails. The extensions of the chart change maps are 
homeomorphisms but not differentiable: This is due to the fact that the Sobolev 
completions Difl['*(M) of Diff(M) of order s > dimM/2 are smooth manifolds 
themselves, but only topological groups: right translations are still smooth, left 
translations and inversions however, are only continuous (and not even Lipschitz). 
So the action of Diff^(M) on Emb*(Af, A^), after Sobolev completion, has aspects 
which are only continuous and thus B''{M,N) = Emb''(M, A^)/Diff'*(M) is only a 
topological manifold in general. This phenomenon also appears in the chart changes 
of the canonical atlas of B{M,N); see [6, 44.1] for an explicit formula of the chart 
change and the role of inversion in Diff''(Af) in it. 



2.5. Covariant curvature and O'Neill's formula, finite dimensional. Let 

p : {E,gE) {B,gB) be a Riemannian submersion between finite dimensional 
manifolds, i.e., for each b € B and x € Ei, := p~^{b) the g^j-orthogonal split- 
ting = T^{Ep(^^)) ® T^iEpf^^))^ =: T^{Epf^^-{) © Hor^(p) has the property that 
TxP : (Hot x{p) T g e) [ThB^gs) is an isometry. Each vector field X G X{E) is de- 
composed as X — -t- X™'' into horizontal and vertical parts. Each vector field 
^ G X{B) can be uniquely lifted to a smooth horizontal field G r(Hor(p)) C 
X{E). O'Neill's formula says that for any two horizontal vector fields X,Y on M 
and any x G E, the sectional curvatures of E and B are related by: 

5p(x)(i?''(p*(^.),P*(>"x))p*(r,),P*(^x)) = g.(i?^(A,,y,)y,,x,)+|||[A,y]--'-|lx- 

Comparing Mario's formula on E and B gives an immediate proof of this fact. Start 
with: 

Lemma. If a £ n^{B) is a 1-form on B, then the vector field {p*af is horizontal 
and we have Tpo [p*a)'^ — a^op. Therefore {p*af' equals the horizontal lift {a^)^°^ . 
For each x d E the mapping (T^p)* : {T*,^B,g^^) {T*E,g^^) is an isometry. 
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Proof. All this holds because for Xx G TxE we have: 

9e{{p*ol)1.,Xx) = {p*a)x{Xx) = ap(^x){Txp.Xx) = ap(^^){Txp.X^°') 

^gE{{p*a)lx'n 
gB{Txp{p*a)i,T,p.Xx) ^ gE{{p*a)lX^°'^) = ap(,)(r,p.X,) 

= 9B{al^^yTxp.Xx). □ 

More generally we have: 

9~E^{p*a,p*l3) = gE{(p*oif, ijp* Pf) = 9B{a^,P'^) o p = p* g^^ {a, 
Consequently, we get for 1-forms a, f3 on B: 

IJE ys 

(p*/3)«|b*af , = )((««)'-) =p*(/3«||af ) 

In the following computation we use 

||[(p*a)«,(p*/3)T°lL=?'*ll[""'/3"]llL 
We take Mario's formula (I2.2p and apply it to the closed 1-forms p*a, p*/3 on E 
where a, (3 are closed 1-forms on B. Using the results above we get: 

45s(i?((p*a)«, b*/3)«)b*/?)«, ip*a)i) = 

= ||d(5^i(p*a,p*/3))||2^. -.g^i(d(|b*a|1^^0,d(|b*/3I|^^0) 

- 3|| [{p*a)K ip*l3)T''\\l - 3||[(/«)«, (P*/3)«]^"||L 

-f 2{p*a)\p*a)H\\p*P\\l-.) + 2(p*/3)«(p*/3)«(|b*a||2 

- 2((p*a)«(p*/3)« + (p*/3)«(p*a)«)g^\p*a,P*/3) 

= p1Mtei(a,/3))f . -p^g5l(d(||af_0,d(||/3f-0) -3p*||[a»,/3«]||^^ 

iJB iiB ^B 

-3\\[ip*a)KipyfV"\\l^ 

+ 2p*(a«a»(||/3||2^0)+2p*(/3»/3«(ll"ll'^0)"2p*((a«/3«+/3V)5B'(«,/3)) 
= 4p*.gs(i?^(a«,/3«)/3«,a«)-3||[(p*a)«,(p*/3)«]™-l|2^ 
which is a short proof of O'Neill's formula. 

2.6. Covariant curvature and O'Neill's formula. Let p : {E,gE) {B,gB) 
be a Riemann submersion between infinite dimensional robust Riemann manifolds; 
i.e., for each b € B and x € E^ := p~^{b) the tangent mapping T^p : {T^E^gE) 
{TbB,gB) is a surjective metric quotient map so that 

(1) \Mgs ■■= inf{^x e TxE : Txp.Xx = 6}- 

The infinimum need not be attained in TxE but will be in the completion T^E. 
The orthogonal subspace {Yx : gE{Yx,Tx{Eb)) = 0} has therefore to be taken in 
T^. 

If Qffc = 5B(a^, ) G gB{TbB) C TJJ'-B is an element in the gs-smooth dual, then 

p*ab := {Txp)*[ab) = gB{al,TxP ) : TxE ^ M is in T*M but in general it is not 
an element in the smooth dual gE{TxE). It is, however, an element of the Hilbert 
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space completion gsiTxE) of the g^-smooth dual gsiTxE) with respect to the norm 
II llg-ii ^-iid the element g]^^{p*ab) =: {p*ai,)^ is in the || ||g^ -completion ExE of 

TxE. We can call g^^(p*ab) —: {p*aby the horizontal lift of = 5i3^(ab) £ TbB. 

In the following we discuss the manifold E and we write g instead of gs- The 
metric gx can be evaluated at elements in the completion TxE. Moreover, for any 
smooth sections X,Y & T{TE) the mapping g{X,Y) : M — R is stiU smooth: 
This is a local question, so let E be c°°-open in a convenient vector space Ve- 
Since the evaluations on X ®Y form a set of bounded linear functionals on the 
space Lgyjjj(yA/; R) of bounded symmetric bilinear forms on Vm which recognizes 
bounded subsets, it follows that g is smooth as a mapping M Lgyjjj(VA/; R), by 
the smooth uniform boundedness theorem, see [6]. 

Lemma. If a is a smooth 1-form on an open subset U of B with values in the 
gB-smooth dual gsiTB), thenp*a is a smooth 1-form on p^^{U) C E with values 
in the \\ W^-i -completion of the gE-smooth dual gE{TE). Thus also {p*a)'^ is 
smooth from E into the gE-completion ofTE, and it has values in the gE-orthogonal 
suhhundle to the vertical bundle in the gE-completion. We may continuously extend 
TxP to the II Wg-i- completion, and then we have Tpo (p*a)'^ = a'^op. Moreover, the 

Lie bracket of two such forms, [{p* a)^ , {p* /3)^ , is defined. The exterior derivative 
d{p*a) is defined and is applicable to vector fields with values in the completion like 
{p*pf. 

That the Lie bracket is defined, is also a non-trivial statement: We have to 
differentiate in directions which are not tangent to the manifold. 

Proof of the lemma. This is a local question; so we may assume that U = B 
and p~^{U) = E are c°°-open subsets in convenient vector spaces Vb and Ve, 
respectively, so that all tangent bundles are trivial. By definition, a' = gB^° 
a : B ^ B X Vb is smooth. We have to show that {p*a)^ = g^^ o p*a is a 
smooth mapping from E into the || ||gg -completion of Ve- By the smooth uniform 
boundedness theorem (see [6]) it suffices to check that the composition with each 
bounded linear functional in a set 5 C is smooth, where 5 C is a set of 
linear functionals on Ve which recognizes bounded subset of Ve ■ For this property, 
functionals of the form gsiv, ) for v G Ve suffice. But 

X i-> {gE)x{v, {p*af\x) = p*a\x{v) = a\x{TxP.v) 

is obviously smooth. 

We may continuously extend the metric quotient mapping TxP to the || ||g^- 
completion and get a mapping TxP : TxE — > TbB where b — p{x). For a second 
form (3 G T{gB{TB)) we have then 

gB{P%,TxP.{p*a)%) - (fib{TxP.{p*a)%) = {p* P)\x{{p* a)%) = 

= gE^{{TxP)*P, {Txp)*a) = gB{Pb,ab) = gB{P^\b, («" o p)[x)) 

which implies Tp o [p*af — o p. 
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For the Lie bracket of two such forms, [(p*a)'', (p* /?)"], we can again assume that 
all bundles are trivial. Then 

[{p*a)K {p*0]{x) = d((p*/3)«)(a;)((p*a)«) - d((p*«)«)(a;)((p*/3)«) 

d((p*/3)«)(x)((p*a)«) = d(.9^i o (Tp)* o /? o p){x){{p*a)») 

= d{g^' o (Tp)* o /3)(6).T,p.(p*a)« 

= d{g^'o (Tp)* o /3){b).aHp{x}). 

So the Lie bracket is well defined. □ 

By assumption, the metric g — Qe admits gradients with respect to itself as in 
(|2.1I 2). in a local chart we have 

(2) D,.zg.{X, Y) = grad^ g{x){X, Y)) = g.{grs.d^ g{x){Z, X),Y) 

for X^Y^Z G Ve- We can then take X.Y ^ Ve in the leftmost expression of (2) 
and thus also in the other two. Then the middle term allows to take Z G VeM 
also. 

Thus the local expressions of the Christoffel symbols of the Levi-Civita covariant 
derivative extend to sections of the completed tensor bundle TE, and therefore 
the Levi-Civita covariant derivative extends to smooth sections of TE which are 
differentiable in directions in TE like {p*a)K Thus expressions like V^.^-jj (p*/?)' 
make sense and are again of the same type so that one can iterate. Thus the 
curvature expression gE{R{{p*oi)^ , {p* (3y){p*a)'^ , {p* f^)"^) makes sense. Moreover, 
all operations used in the proof of p.2p work again, so this result holds. The proof 
in (|2.6p works and we can conclude the following result: 

Theorem. Let p : {E,gE) — > {B,gB) be a Riemann submersion between infinite 
dimensional robust Riemann manifolds. Then for 1-forms a, {3 ilg^{B) O'Neill's 
formula holds in the form: 

5B(i?^(a«,/3«)/3«,a«) 

= gE{R'^{{p*a)\ {p*l3)^){p*l3)K + {P*(3)TX, 



3. The diffeomorphism group Diff5(Af) 

3.1. DifFeomorphism groups. Let TV be one of the following: 

• is a compact manifold: Then let DifF(Af) be the regular Lie group [6l 
section 38] consisting of all smooth diffeomorphisms of M. 

• N is M": we let Diff5(R") denote the group of all diffeomorphisms of R" 
which decay rapidly towards the identity. This is a regular Lie group (for 
n = 1 this is proved in [TOl 6.4]; the proof there works for arbitrary n). 
Its Lie algebra is the space Xs(]R") of rapidly falling vector fields, with the 
negative of the usual bracket as Lie bracket. 

• More generally, (A'^, g) is a non-compact Riemannian manifold of bounded 
geometry, see [3]. It is a complete Riemannian manifold and all covariant 
derivatives of the curvature are bounded with respect to g. Then there is a 
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well developed theory of Sobolev spaces on N; let H°° denote the intersec- 
tion of all Sobolev spaces which consists of smooth functions (or sections). 
Even on = M the space H°° is strictly larger than the subspace S of 
all rapidly decreasing functions (or sections) which can be defined by the 
condition that the Riemannian norm of all iterated covariant derivatives 
decreases faster than the inverse of any power of the Riemannian distance. 
There is nearly no information available on the space S for a general Rie- 
mannian manifold of bounded geometry. For the following we let S denote 
either H°° or the space of rapidly decreasing functions. We let DifF5(A^) 
denote the group of all diffeomorphisms which decay rapidly towards the 
identity (or differ from the identity by H°°). It is a regular Lie group with 
Lie algebra the space Xs{N) of rapidly decreasing vector fields with the 
negative of the usual bracket. In [10, 6.4] this was proved for N = M., but 
the same proof works for the general case discussed here. 

In general, we need to impose some boundary conditions near infinity for groups 
of diffeomorphisms on a non-compact manifold M: The full group Diff(A^) of all 
diffeomorphisms with its natural compact C°° topology is not locally contractible, 
so it does not admit any atlas of open charts. 

For uniformity of notation, we shall denote by Diff5(A^) any of these regular 
Lie groups. Its Lie algebra is denoted by Xs{N) in each of these cases, with the 
negative of the usual bracket as Lie bracket. We also shall denote by O = C°° fl S' 
the space of smooth functions in the dual space S' (to be specific, this is the space 
Oi\i in the sense of Laurent Schwartz, if = R"). 

3.2. Riemann metrics on the diffeomorphism group. Motivated by the 
concept of robust Riemannian manifolds and by [16] chapter 12] we will construct a 
right invariant weak Riemannian metric by assuming that we have a Hilbert space 
Ti together with two bounded injective linear mappings 



where T(j2{TN) is the Banach space of all vector fields X on N which are 
globally bounded together with V^X and V^V^X with respect to g, such that 
j2 o ji '■ (TN) — > r(j2 (TN) is the canonical embedding. We also assume that ji 
has dense image. 

Dualizing the Banach spaces in equation (1) and using the canonical isomor- 
phisms between T-L and its dual Ti' - which we call L and we get the diagram: 



(1) 



Xs{N)=rs{TN) 



> n {TN) 



(2) 



Ts{TN) 



Ts'{T*N) 



31 



L 



K 



Tm2{T*N) 
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Here we have written rsi{T*N) for the dual of the space of smooth vector fields 
Ts{TN) = Xs{N). We call these 1-co-currents as 1-currents are elements in the 
dual of rs{T*N). It contains smooth measure valued cotangent vectors on N 
(which we will write as rs{T*N vol(A^))) and as well as the bigger subspace of 
second derivatives of finite measure valued 1-forms on N which we have written as 
r^/s {T*N) and which is part of the dual of (TN). In what follows, we will have 
many momentum variables with values in these spaces. 

The restriction of L to Xs{N) via ji gives us a positive definite weak inner 
product on Xs{N) which may be defined by a distribution valued kernel - which 
we also write as L: 

{ , )l:Xs{N) X XsiN) -J> R, defined by 

J JNxN 

where L £ rs'(prt(T*iV) pr*(r*A^)) 

Extending this weak inner product right invariantly over Diff5(A'^), we get a robust 
weak Riemannian manifold in the sense of 12.41 

In the case (called the standard case below) that N = R" and that 
{X,Y)l= I {{l-AAyx,Y)dx 

we have 

^(^'2/) = (7^ / e'^^^^-yHl + m'ydAj2{du%^dx)^idu\®dy) 
\ (27r)" 7^gR„ / 

where d^, dx and dy denote Lebesque measure, and where (u*) are linear coordinates 
on R". Here Ti. is the space of Sobolev i?' vector fields on N. 

Note that given an operator L with appropriate properties we can reconstruct 
the Hilbert space H with the two bounded injective mappings ji , j2 . 

Construction of the reproducing kernel K : The inverse map K is even nicer as it 
is given by a tensor, the reproducing kernel. To see this, note that T f,.[i{T*N) 
contains the measures supported at one point x defined by an element Oix G T*N . 
Then i2{K{i'2{ax))) is given by a vector field K^^ on A^which satisfies: 

(3) {K^^,X)H=a.,{32X){x) for all X e H, a, £ T^N. 

The map ax > is weakly C^, thus by ^6, theorem 12.8] this mapping is strongly 
Lip^ (i.e., differentiable and the derivative is locally Lipschitz, for the norm on H). 
Since eVyoK : T*N B ax Ka^{y) G TyN is linear we get a corresponding 
element K{x,y) € L(T*N,TyN) = TxN (S)TyN with K{y,x){ax) = Ko,^{y). 

Using (3) twice we have (omitting J2) 

Py.K{y,x){ax) = {K{ ,x){ax),K{ ,y){l3y))-H ^ ax-K{x,y){/3y) 

so that: 

. Kix, yY = K{y, x) : T^N ^ T.TV, 
• K e Tci (pri* TN ® pr^ TN). 
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Moreover the operator defined directly by K by integration 



K : Tm^{T*N) ^Ts{TN) 



JyieN 



is the same as the inverse K to L. In fact, by definition, they agree on sections 
in r(72 (T*M) with finite support and these are weakly dense. Hence they agree 
everywhere. 

We will sometimes use the abbreviations {a\K\, \K\I3) and {a\K\l3) for the con- 
traction of the vector values of K in its first and second variable against 1-forms a 
and 0. Often these are measure valued 1-forms so after contracting, there remains 

a measure in that variable which can be integrated. 

Thus the tensor K determines L and hence H and hence the whole metric 
on DiSs{N). It is tempting to start with the tensor K, assuming it is symmetric 
and positive definite in a suitable sense. But rather subtle conditions on K are 
required in order that its inverse L is defined on all infinitely differentiable vector 
fields. For example, if N = W, the Gaussian kernel K{x,y) = e"!^""! does not 
give such an L. 

In the standard case we have 



where Ki is given by a classical Bessel function of differentiability class C^' . 

3.3. The zero compressibility limit. Although the family of metrics above 
do not include the case originally studied by Arnold - the metric on volume 
preserving diffeomorphisms - they do include metrics which have this case as a 
limit. Taking N = W^ and starting with the standard Sobolev metric, we can add 
a divergence term with a coefficient B: 



Note that as B approaches oo, the geodesies will tend to lie on the cosets with 
respect to tlic subgroup of volume preserving diffeomorphisms. And when, in ad- 
dition, A approaches zero, we get the simple metric used by Arnold. This 
suggests that, as in the so-called 'zero- viscosity limit', we should be able to con- 
struct geodesies in Arnold's metric, i.e. solutions of Euler's equation, as limits of 
geodesies for this larger family of metrics on the full group. 

The resulting kernels L and K are no longer diagonal. To L, we must add 
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It can be checked that the corresponding kernel K will have the form 

"99"" d d 

K{x, y) = K,{x-y)Y, — ® —+Y, -y^d^^®d^ 

i=l ^ i=l j=l ^ 

where Kq is the kernel as above for the standard norm of order I and Kb is a second 
radially symmetric kernel on M" depending on B. 



3.4. The geodesic equation. According to [2], the geodesic equation on any 
Lie group G with a right-invariant metric is given as follows. Let g(t) be a path in 
G and let u{t) — g{t).g{t)~^ — r(/i^*^*-' )g{t) be the right logarithmic derivative, a 
path in its Lie algebra g. Here fj,^ : G —?' G is right translation by g. Then g{t) is a 
geodesic if and only if 

dtu = — ad^ u. 

where the transposed adj^ is the adjoint of a,dx '■ Q ^ B with respect to the metric 
on 0. 

In our case the Lie algebra of Diff 5 (N) is the space X5 (N) of all rapidly decreas- 
ing smooth vector fields with Lie bracket (we write adx Y) the negative of the usual 
Lie bracket a.dx Y = — [X,Y]). Then a smooth curve t (p{t) of diffeomorphisms 
is a geodesic for the right invariant weak Riemannian metric on Diff5(Af) induced 
by the weak inner product ( , )l on Xs{N) if and only if 



dtU 



ad' 



as above. Here the time dependent vector field u is now given by dt^{t) — u{t)o (p{t), 
and the transposed adj is given by 

(adj r, Z)l = (y, adx Z)l - -{Y., [X, Z])l. 

The inner product is weak; existence of adj is equivalent to condition (1) for ro- 
bustness of the weak Riemannian manifold (Diff 5 (iV), ( , This is equivalent 
to the fact that the dual mapping ad^ : 'X.s{N)' — >■ Xs^N)' maps the smooth dual 
L{Xs{N)) to itself; we also have L o adj = a.d*xoL. Using Lie derivatives, the 
computation of ad^ is especially simple. Namely, for any section u of T*N (g) vol 
and vector fields G Xs{N), we have: 



IN Jn Jn 

hence ad^(w) = +C^(uj). Thus the Hamiltonian version of the geodesic equation 
on the smooth dual L{Xs{N)) C rc2{T*N (g) vol) becomes 



dta 



■ ad 



K(a)< 



= -C 



or, keeping track of everything. 



(1) 




One can also derive the geodesic equation from the conserved momentum mapping 
J : TBiSsiN) s' XsiN)' given by Jig,X) = Lo KdigY X where kd{g)X = Tgo 
X o g^^. This means that Ad{g{t))u{t) is conserved and — dt Ad{g{t))u{t) leads 
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quickly to the geodesic equation. It is remarkable that the momentum mapping 
exists if and only if (Diff5(A^), ( , )l is a robust weak Riemannian manifold. 



4. The differentiable Chow manifold (alias the non-linear 

Grassmannian) 



4.1. The differentiable Chow manifold as a homogeneous space for 
DifF5(A^) and the induced weak Riemannian metric. Let M be a compact 
manifold with dim(A/) < dim(A^). The space of submanifolds of N diffeomorphic 
to M will be called B{M, N). In the case m = and iV = R^, i.e. M is a finite set 
of, say p, points in Euclidean D-space, the space B{M,N) is what we called the 
space of landmark points C{M.^) in our earlier paper [8]. 

B{M, N) can be viewed as a quotient of DifF5(iV). If we fix a base submanifold 
Fq C N diffeomorphic to M, then we get an action of Diff5(iV) on B{M,N) by 
(fi I—)- (p{Fo). The image will be an open subset Bo{M, N) of B{M, N) which is the 
quotient of Diff5(A^) by the subgroup of diffeomorphisms which map Fq to itself. 
We will study B{M, N) using this approach and without further comment replace 
the full space B{M, N) by this component Bo{M, N). 

The normal bundle to F C iV may be defined as TB^ C TN\b, with the help 
of an auxiliary Riemann metric on N. But we want to avoid this auxiliary metric, 
so we shall define the normal bundle as the quotient Nor(F) := TN\f/TF over 
F. Then its dual bundle, the conormal bundle, is Nor*(F) = Annihilator(TF) C 
T*N\f, a sub-bundle not a quotient. The tangent space TpB^M, N)io B{M, N) at 
F can be identified with the space of all smooth sections r5(Nor(F)) of the normal 
bundle. 

A simple way to construct local coordinates on B{M,N) near a point F G 
B{M, N) is to trivialize a neighborhood oi F C N. To be precise, assume we have 
a tubular neighborhood, i.e., an isomorphism $: 



B{M,N) Nor(F) 
U U 

Ub Un 

u u 

F ~ 0-section 



from an open neighborhood Ub of F in N to an open neighborhood Un of the 
0-section in the normal bundle Nor(i^). Assume moreover that $ is the identity on 
F and its normal derivative along F induces the identity map on Nor(F). The map 
$ induces a local projection n : Ub ^ F and partial linear structure in the fibres of 
this projection. Then we get an open set J7$ C B{M, N) consisting of submanifolds 
F' C Ub which intersect the fibres of tt normally in exactly one point. Under $ 
these submanifolds are all given by smooth sections of Nor(F) which lie in Un- If 
we call this set of sections Ur we have a chart: 



B{M,N) DU^^UrC rs(Nor(F)) 
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We define a Riemannian metric on B{M,N) following the procedure used for 
the space of embeddings. For any F c N, we decompose H into: 

•T^vert = j2^{{X e Tc^iTN) : X{x) G T^F, for aU x G F}) 

H^'^ = perpendicular complement of 'H^'^^ 

It is then easy to check that we get the diagram: 



Ts{TN)^ 



TsCNoriF))^ 



^horc 



rc=(riv) 



rc.(Nor(F)). 



As this is an orthogonal decomposition, L and K take 7^™"^* and 7^^°"" into their 
own duals and, as before we get: 



r5(Nor(F)) 



r5'(Nor*(F)) 



31 



ftp 



Kf 



32 



Tc2{Nor{F)) 



rM2(Nor*(F)) 



Kp is just the restriction K to this subspace of H' and is given by the kernel: 

Kp{xi,X2) := image of K{xi,X2) € NoTxi{F) (8> Nor^^^ (F)), xi,X2 € F. 

This is a section over F x F of pr* Nor(F) (g) pr2Nor(F). We can identify 
the space of horizontal vector fields 7^^°"^ as the closure of the image under Kp 
of measure valued 1-forms supported by F and with values in Nor*(F). A dense 
set of elements in H^'^ is given by either taking the 1-forms with finite support or 
taking smooth 1-forms. In the first approach, TL^'^ is the closure of the span of the 
vector fields |i4ri?(-,x)|aa:) where x G F and £ Nor*(F). In the smooth case, fix 
a volume form k on M and a smooth covector ^ G r5(Nor*(F)). Then ^.k defines 
a horizontal vector field h like this: 



h{xi) 



= / \Kp{xi,X2)\^{X2).k{X2)) 

Jx2eF 



The horizontal hft h^"'' of any h G TpB{M, N) is then: 



1-2 GF 



\K{yi,X2)\LFh{x2)), 1/1 e TV. 



Note that all elements of the cotangent space a G r5/(Nor*(F)) can be pushed up 
to N by (j>)*, where jp : F ^ N is the inclusion, and this identifies {jp),,a with 
a 1-co-current on N. 

Finally the induced homogeneous weak Riemannian metric on B{M, N) is given 
like this: 

{h,k)p= f (/i^-(yi),L(A:^-)(2/i))= / {K{Lph)){y,),{Lpk){y,)) 
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-L 



(if (2/1, 2/2), {LFh){yi) ® {LFk){y2)) 



(!/l,!/2)eJVxiV 

{LFh{x-i_)\KF{xi,X2)\LFh{X2)) 

^ {xi,xi)eFxF 

With this metric, the projection from Diff5(A^) to B{M, N) is a submersion. The in- 
verse co-metric on the smooth cotangent bundle \_\f^b{m n) r(Nor*(i^)(8)vol(F)) c 
T*B{M, N) is much simpler and easier to handle: 



{a,(3}F 



FxF 



{a{xi)\KF{xi, X2)\a{xi)) . 



It is simply the restriction to the co-metric on the Hilbert sub-bundle of T* DifF5(iV) 
defined by n' to the Hilbert sub-bundle of subspace T*B{M, N) defined by Tip. 

Because they are related by a submersion, the geodesies on B{M,N) are the 
horizontal geodesies on DifF5(A^), as described in box (1), section 3.4. We have two 
variables: a family {Ft} of submanifolds in B{M, N) and a time varying momentum 
a{t, •) G Nor*(Fi) eg) vol(Ft) which lifts to the horizontal 1-co-current {jFt)*{oi{t, ■) 
on N. Then the horizontal geodesic on DifF5(A^) is given by the same equations as 
before: 



dt{Ft) 


= resNor(Fo(w(*r)) 


u{t,x) 


= / \K{x,y)\a{t,y)) eXs{N) 






dtiijFMc^it,.)) 


= -^u{t,-)iijFt)*i(^it,-))- 



This is a complete description for geodesies on B{M, N) but it is not very clear how 
to compute the Lie derivative of {jFt)*{oi{t, ■). One can unwind this Lie derivative 
via a torsion-free connection, but we turn to a different approach which will be 
essential for working out the curvature of B{M, N). 



4.2. Auxiliary tensors on B{M,N). Our goal is to reduce calculations on the 
infinite dimensional space B{M,N) to calculations on the finite dimensional space 
N. To do this we need to construct a number of useful tensors on B{M, N) from 
tensors on TV and compute the standard operations on them. These will enable us 
to get control of the geometry of B{M, N). Let m be the dimension of M, n the 
dimension of N. For F e B{M,N), let Jf ■ F ^ N he the embedding. We will 
assume that M is orientable for simplicity, so that vol(M) = 0'"(M). 

(1) We denote by i the left action: 

£ : Diff5(A/') X B{M, N) B{M, N) 

given by £{(p,F) or i^{ip) = <fi{F). For a vector field X e Xs{N) let Bx be the 
infinitesimal action (or fundamental vector field) on B{M,N) given by Bx{F) = 
Tid{i^)X with its flow Flf''(F) = Flf (F). The fundamental vector field mapping 
of a left action is a Lie algebra anti-homomorphism and the Lie bracket on Diff5(A'') 

is the negative of the usual Lie bracket on Xs{N), so we have [Bx, By] = B[x,y]- 
The set of these vectors {Bx{F) : X e Xs{N)} equals the whole tangent space 
TfB{M, N). 
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(2) Note that B{M, N) is naturally submanifold of the vector space of m-currents 
on N: 

B{M, N) Ts'{n"'{N)), yia.F^ (^Lj^ J u?j . 

Any a G fl"^{N) is a linear coordinate on Ts'{^"^{N)) and this restricts to the 
function B^ S C°°{B{M,N),R) given by = a. If a = d/S for /3 e 

n"^-\N) then 

= Ba0{F) = [ fpdp = [ dfpp = 
Jf J f 

by Stokes' theorem. 

For a G Qr^{N) and X e Xs{N) we can evaluate the vector field -B^ on the 
function B^- 

Bx{B^){F) = dB^{Bx){F) = d,\oB„{Fl^ {F)) = [ f^Cxa = Bc,(^){F) 

Jf 

as well as = / jp{ixda + dixa) = / jpixda = 
Jf Jf 

If X e Xs{N) is tangent to F along F then Bx{Ba){F) = Cx\p3*pa = 0. 

More generally, a pm-form a on A''*^ defines a function B^'' on S(M, N) by 
Ba\F) = a. Using this for p = 2, wc find that for any two m-forms a, p on 
N, the inner product of B^ and -B^ is given by: 

Qg^iB^B^) = -B[^|k|^)- 

(3) For a G 12™+''(iV) we denote by B^ the fc-form in n''{B{M,N)) given by the 
skew-symmetric multi-linear form: 

{B,,)FiBxAF), ■ ■ ■,BxjF)) = f jF*iix,A-AX,a). 

Jf 

This is well defined: If one of the Xi is tangential to at a point x G F then jp* 
pulls back the resulting m-form to at x. 

Note that any smooth cotangent vector a to F G B{M,N) is equal to Ba{F) 
for some closed [m + l)-form a. Smooth cotangent vectors at F arc elements of 
rs(F,Nor*(F) (g) n'^{F)). Fix a nowhere zero global section k of f^™(F). Then ^ 
is the differential of a unique function / on the normal bundle to F which is linear 
on each fibre. Let be a local isomorphism from a neighborhood of F in A/^ to a 
neighborhood of the 0-section in this normal bundle and let p be a function on the 
normal bundle which is one near the 0-section and has support in this neighborhood. 
Take a = d{f.K o (p) (extended by zero). It's easy to see that this does it. 

Likewise, a pm + k form a G 17^""+*^ (A^^) defines a fc-form on B{M, N) as follows: 
First, for X G Xs{N) let g X{Np) be given by 

^inl,...,np) ■= (^"1 X On^ X • • • X 0„J + (0„, X X 0„3 X . . . X 0„J + . . . 

■•■ + (0„i X ... X 0„^_, X X„J. 
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Then we put 



This is just B applied to the submanifold F'p C and to the special vector fields 
Thus all properties of B continue to hold for iJ^^^; in particular, (4) below 
hold for instead of X. 



(4) We have 



because 



(«Bxi Ba) {Bx2 , • • • , Bx^ ) (F) = Ba , Bx2 Bx^ ) if) 

= j^jF*{ixk---ix2{ix,cx)) = Bi^^a{Bx2,---,Bx^){F) 

For the exterior derivative we have 



dBc, = Bdoc for any a G f2™+'=(A/'). Namely, 



{dB^){Bxo,-,Bx,KF) = J2i-^yBxABa{Bxo,-,B^,,-,BxMF) 

i=0 



%<3 



j]AXo/\--Xi--Xj--/\Xk 



= E(-ir/ j*Fix,dixoA-x.-r.x,^ + ll^-^y^' I fF\x.,x,] 

i=0 ■'F i<j ''F 

^ i=0 

~ E^"-*-) ^XoA-%-AXj-^A[Xi,Xj]AXj+-i_-AXu)^ 
i<j 

= / jF*^{-lT(^CxJx^- -ix,- -ixo 

•'F 1=0 

k 

•'F 4=0 

= / • ■^Xi+^dix, ■ -ixo + ^{-'^yixk- -ix^dix^.^- -1x0)0: 

•'F i=0 i=0 

= + y" jF*iXk ■ -ixoda = Bda{Bxo,--,BxJ{F) 



a 



Finally we have CBxBa = Bc^ 



since 



^BxBa — i^Bx d + diBx)Ba — -B(ix(i+dix)a — Bcxc 
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Note that these identities gcnerahzc the resuhs in item (2). 

(5) For a e Vf^+'^{N) we pull back to Diff5(7V) the 1-form on B{M,N) where 
^o(fb) = F: 



{{e''-YB^)jXo^,) = {{l^YB^)^^{X) = {B^)f{Bx{F)) = / rp^x 



F 



fi : n"'+\N) X B{M,N) XsiN)' 

^{a, F) is a 1-cocurrcnt with support along F. 



The mapping ^ : Qr+^{N) x B{M,N) XdN)' is smooth, : 17"+i(A^) ^ 
Xc(iV)' is bounded linear, and the differential of /j-a ■ B{M,N) 3Cs{Ny is com- 
puted as follows: 

{d{fia){Bx{F)),Y) = {DF,B^f,{a,F),Y) = DF,B^{fx{a,F),Y) = dt\o{aFix^p^,Y) 



= dt\o [ JFif (F)*^Ya = dt\o f (Plfojp o (Flf \F)-^riYa 
= dt\o[ {Flf If)-')* jF*{FlfriYa 

Jpif (F) 

= dt\o / jF*(Flf )*iYa = / jp* jCxiiya) 

J F J F 

= J^jF*{i[x,Y]a + iYjCxa) = {fx{a,F),jCxY) + {fx{jCxa,F),Y). 



This means 



(6) dfi„{Bx(r)) = fi[n.F) o Cx +f,{Cx<^.F) = -C.xfi{o. F) + fi{C.xo. F), 



where Cxl^{oi, F) denotes the Lie derivative of 1-currents. There are two interpre- 
tations of formula (6): 

dfJ.a{Bx) = -Cx O HCxa, 

dtio,{Bx{F)) = -[Cx^i^){a). 

We shall also need the mapping : ^'^{N) x B{M,N) C^{N)' with values in 
the linear space of distributions (without the density part) on N which is given by 



(/x(7,i^),/)= / /•7= / 3F*{91). 

JF J F 



The distribution iJ,{'y,F) is again bounded linear in 7 S ri'"(7V), and its derivative 
with respect to F is given by (6) again, with the same proof as above. 

5. GeODESICS AND CURVATURE ON B{M,N) 

We want to use the auxiliary tensors of the last section to derive formulas for 
geodesies and curvature on B{M), using Mario's formula to compute the curvature. 
The basic idea is to write a smooth co- vector a at a point F G B{M,N) as Ba 
where a is an (m + l)-form on N. As always, for any (m + l)-form a on N, B^ is 
the (C^) vector field on B{M, N) which is dual to the smooth 1-form Ba- At each 
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point F E B, lifts horizontally to a tangent vector at the identity to Diff5(iV), 
which is given by the vector field 

Jn 

so that = Bi{F). See K^5). 

With these co-vectors, we consider next the force introduced in section 2.3. We 
have: 

2^^(a,/3) 

But {a\K\f3) is a 2m-forin on N x N and d can be split into two parts di+d2 acting 
on the first and second factors. Evaluating this 1-form at F and taking its inner 
product with Bx,X G Xs{N), we get: 

= 11 jFxF*i{ix)i{di{{a\K\m + {^x)2{d2{{a\K\m) 
J Jfxf 

because F x F has type (to, to) and the integrand must have the same type 
jpixd (v(0,F)«(a) + 
hence 

2F{a, 13) = = 5^, 7 = '^^m(/3,f)« («) + ^f.{a,F)» (/?)■ 

Here the superscript 2 on X means that X'^^^ is the vector field on N x N given by 
X X + X X whereas on B, because d{{a\K\P)) is a (2to + l)-form on x iV, 
we must apply B^'^\ not B, to it. Thus we define the force F using operations on 
the finite dimensional manifold N by: 



FN{a,P,F) := (image in Nor*(i^)«.vol(F)) (^^(/:,.(^,f (a) + (/3)) j . 

The term 'force' comes from the fact that the geodesic acceleration is given by 
F{a,a). In our case, we find that the geodesic equation on B{M,N) can be ex- 
tended to an equation in the variables Ft € B{M,N) and a{t,-) a time varying 
(to -|- l)-form on N: 



dt{Ft) 


= (res to Nor(F))u 


u 


= M(a,F)»= / \K{;y)\a{y)) 






dt{a) 


= F{a,a,F) ^ Cuia). 



Moving to curvature, fix F. Then we claim that for any two smooth co- vectors 
a, b at F, we can construct not only two closed (to -|- l)-forms a, /3 on iV as above 
but also two commuting vector fields Xa , Xfj on N in a neighborhood of F such 
that: 

(1) Ba{F) = a and B/3(F) = b, 

(2) BxAF) = a" and Bx,{F) = fe» 
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(3) Hx^ (a) = Lx^ (/3) = Lx, (a) = Hx, (/?) = 

(4) [X„,X^]=0 

We can do this using a local isomorphism of N with the normal bundle to F in iV 
as above. This gives a projection vr of a neighborhood of F in N to F and partial 
linear structure on its fibres. Then for a and /? use (m + l)-forms k A uj where k 
is a pull back of an m-form on F and w is a 1-form constant along the fibres; and 
for Xa and use vector fields which are tangent to the fibres of tt and constant 
with respect to the linear structure on them. 



We are now in a position to use the version (|2.3p of Mario's formula. As it stands, 
this formula calculates curvature using operations on B{M,N). What we want to 
do is to write everything using forms and fields on N instead. We first need an 
expression for the stress D{a,l3) in this formula. Using notation from p. 31 2): 



Dia,P,F)=D 



F,Bx^(F)y-Dp 



= [Bx^ , Bl - Bx, ] (F) = [Bx^ , i?» ] (F) . 

In order to compute the Lie bracket, we apply it to a smooth function B^ on 
B{M,N) where 7 £ 17™(A^). Then we have, using S^l repeatedly: 



,«7 



{C^,B,){F) = {Cb^^^^^^^,B,){F)^Bc^^^ 

{^Bx^^B*^B-l)^F) ^ {Cbx Bl 

= B{CDj.s^^KP.F)n)\.J' ) + oyL.x^i-i,(f}^F)r 
{Cs^CBxMiF) = {Cb^^^,^,Bcx^.,)(F) = B{C,(p^F)^Cx^-i){F) 



n){F) + B{Cx^c^(^p^F)n){F) 



Dfb 



\KHP,F)) 



N 



\K\D 



F,B^ 



AP^F)) 



N 



f \K\{-Cx^^i{P,F)+^iiCx^p,F))) bygjle) 

JN 



\Coxx^K\^iil3,F)) + fi{Cx^p,F)^ 



N 

{[Bxa, Bf,(^,,^ p)t]By){F) = {CgiB^f - Cgt CBx^By){F) 
= Bi£o,,s,^^.i(i,F)n)iF) + Biqx^^^.(p,F)m){F) 

= {^B(Df,Bx„ P^{I3,F)» + [X„MI3,F)»])B',){F) 

[Bx^ , B^ibe^F)t]iF) - BiDF^Bx^KP. F)i + CxMP. Ff) 

^b( f \Co>,x^K\fi{l3,F))+fi{Cx^p,F)^ + + f \Cx^^oK\fiW, F))) 

= B^jjCxj^)KAp,F)))+Q. 



Thus we define the stress on A'' by: 



DN{a,l3,F){x) = (restr. to Nor(F)) 



y&F 
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Next consider the second derivative terms in i?ii. A typical term works out as 
follows: 

Bx„BxA< Bp, Bp >) = Cbx^^Bx^{< Bp, Bp >) = Bc^^^^c^^2)<^3\K\p> 

= B<p\Cn^^^C^^2)K\p> 

Extending Lie bracket notation slightly, we can write 

Analogous formulas old for the other terms. 

Finally, putting everything together, we find the formula for curvature: 



/3 



{RB{M,N){BlBj,)BlBi){F) = Rn + R12 + R2 + R3 



Rii 



FxF 



(^{I3\C^(2)C^(2)K\(3) + (a\C^(2) C^(2) K\a) 



R12 — 



{D{a, a, F),F{(3, (3, F)) + {D{(3, (3, F), F{a, a, F)) 
- {D{a, /3, F) + D{(3, a, F), F(a, /3, F))) 



i?2 = \\F{a,p,F)\\l^ - {F{a,a,F)),F{l3,l3,F)) 
i?3 = -l\\DNia,p,F) - DN{f3,a,F)\\l^ 



Kb 



In the case of landmark points, where m = 0, N = M.^ and K is diagonal, it is easy 
to check that our force and stress and the above formula for curvature are exactly 
the same as those given in our earlier paper [S] . In that paper the individual terms 
are studied in special cases giving some intuition for them. 



6. Appendix on Convenient Calculus - Calculus beyond Banach 

SPACES 

The traditional differential calculus works well for finite dimensional vector 
spaces and for Banach spaces. For more general locally convex spaces we sketch 
here the convenient approach as explained in [5] and [6]. The main difficulty is 
that composition of linear mappings stops being jointly continuous at the level of 
Banach spaces, for any compatible topology. We use the notation of [Fj and this is 
the main reference for the whole appendix. 



6.1. Convenient vector spaces and the c°°-topology. Let i? be a locally con- 
vex vector space. A curve c : M ^> i? is called smooth or C°° if all derivatives 
exist and are continuous - this is a concept without problems. Let C°°(R,£') be 
the space of smooth functions. It can be shown that C°°{M.,E) does not depend 
on the locally convex topology of E, but only on its associated homology (system 
of bounded sets). 
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E is said to be a convenient vector space if one of the following equivalent con- 
ditions is satisfied (called c°°-completeness): 

(1) For any c G C°°(K, E) the (Riemann-) integral c{t)dt exists in E. 

(2) A curve c : R — ?> £' is smooth if and only if A o c is smooth for all A G E', 
where E' is the dual consisting of all continuous linear functionals on E. 

(3) Any Mackey-Cauchy-sequence (i. e. tnm{xn — Xm) for some tnm oo 
in M) converges in E. This is visibly a weak completeness requirement. 

The final topology with respect to all smooth curves is called the c°°-topology on 
E, which then is denoted by c°°E. For Frechet spaces it coincides with the given 
locally convex topology, but on the space V of test functions with compact support 
on M it is strictly finer. 

6.2. Smooth mappings. A curve c : M — > i? is called smooth or C°° if all deriva- 
tives exist and are continuous - this is a concept without problems. Let C°°(K,i?) 
be the space of smooth functions. It can be shown that the set C°°(K, E) does not 
depend on the locally convex topology of E, only on its system of bounded sets. 
Let E, F, and G be convenient vector spaces, and let C/ C E be c°°-open. Here 
is the key definition that makes everything work: a mapping f : U ^ F is called 
smooth or C°°, if / o c e C°°(R, F) for ah c e C°°(M, U). 

The main properties of smooth calculus are the following. 

(1) For mappings on Frechet spaces this notion of smoothness coincides with 
all other reasonable definitions. Even on this is non-trivial. 

(2) Multilinear mappings are smooth if and only if they are bounded. 

(3) If f : E D U ^ F is smooth then the derivative df : U x E F is smooth, 
and also df : U L(E,F) is smooth where L(E,F) denotes the space of 
all bounded linear mappings with the topology of uniform convergence on 
bounded subsets. 

(4) The chain rule holds. 

(5) The space C°°{U,F) is again a convenient vector space where the structure 
is given by the obvious injection 

C--{U,F)^-^^ -Q c°°{R,^), f^{lofoc)c,i, 

where (M, M) carries the topology of compact convergence in each deriv- 
ative separately. 

(6) The exponential law holds: For c°°-open V (Z F , 

C°^{U,C°°{V,G)) ^ C°°{U X V,G) 

is a linear diffeomorphism of convenient vector spaces. Note that this is the 
main assumption of variational calculus where a smooth curve in a space 
of functions is assumed to be just a smooth function in one variable more.. 

(7) A linear mapping f : E ^ C°°{V,G) is smooth (bounded) if and only if 

E ^ > C°°(y,G) > G is smooth for each v G V. This is called the 

smooth uniform boundedness theorem ^ 5.26]. 
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(8) The following canonical mappings are smooth. 

ev : C°° {E, F) x E ^ F, ev(/, x) = f{x) 
ins-.E^ C°°{F,E X F), ins(x)(y) = {x,y) 
( ■.C'^{E,C°°{F,G)) -^C^iE X F,G) 
( : C°°(i; X ^^,6) C°°(i;,C°°(i^,G')) 
comp : C°°{F,G) x C°"{E,F) C°°{E,G) 

C°°( , ):G'^(F,Fi)xG°°{Ei,E)^G'^(G'^{E,F),C°°{Ei,Fi)) 
(/,.9) ^{h^ fo ho g) 

n : n ^"(^- ^ ^"(11 n 

Smooth mappings are always continuous for the c°°-topology but there are 
smooth mappings which are not continuous in the given topology of E. This is 
unavoidable and not so horrible as it might appear at first sight. For example the 
evaluation E x E* ^ R is jointly continuous if and only if E is normable, but it is 
always smooth. 
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